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ABSTRACT 



We perform a global fit to the inclusive structure function considering a QCD inspired model 
based on the summation of gluon ladders describing the ep scattering. In lines of a two Pomeron 
approach, the structure function F2 has a hard piece given by the model and the remaining soft 
contributions: the soft Pomeron and non-singlet content. We have investigated several choices 
for the soft Pomeron and its implication in the data description. In particular, we carefully 
estimated the relative role of the hard and the soft contributions in a large span of x and Q^. 



1 Introduction 



A great challenge in understanding the hard interactions has been posed by the HERA small-x 
data [1] in order to describe the strong growth of the inclusive structure function F2 as the 
Bjorken scale x decreases. This feature is also supplemented by the scaling violations on the 
hard scale given by the photon virtuality Q^. The currently high accuracy reached turns out the 
measurements having quite small statistic uncertainties, probing a large interval of virtualities, 
i.e. from a real photon to thousands of GeVs. As far the Regge approach is concerned, at high 
energy the ep scattering process is dominated by the exchange of the Pomeron trajectory in the 
t-channel. From the hadronic phenomenology, this implies that the structure function would 
present a mild increasing on energy (s ~ Q^/x) since the soft Pomeron intercept ranges around 
ajp(O) ~ 1.08. Such behavior is in contrast with high energy ep data, where the effective 
intercept takes values Aeg — 1-3 — 1.4. In the Regge language this situation can be solved 
by introducing the idea of new poles in the complex angular momentum plane, for instance 
rendered in the multipoles models [2, 3, 4], producing quite successful data description. Other 
proposition is given by the two Pomeron model [5], introducing an additional hard intercept and 
corresponding residue. However, a shortcoming from these approaches is the poor knowledge 
about the behavior on virtuality, in general modeled in an empirical way through the vertex 



functions. An interesting mixed procedure is to fix the initial conditions for the QCD parton 
distributions (valence and sea/gluon) in a sufficiently large initial scale Qq^1 — 5 GeV^ from 
the Regge phenomenology and then perform the QCD evolution up to higher virtualities [6, 7]. 

On the other hand, the high photon virtuality allows the applicability of the QCD pertur- 
bative methods. Two main approaches emerge from the QCD formalism: the DGLAP and 
BFKL evolution schemes. The DGLAP formahsm [8] is quite successful in describing most of 
the measurements on structure functions at HERA and hard processes in the hadronic collid- 
ers. This feature is even intriguing, since its theoretical limitations at high energy arc well 
known [9]. Other pcrturbative approach is the BFKL formalism [10], well established at LO 
level but not yet completely understood at NLO accuracy. The main issue in the NLO BFKL 
effects is the correct account of the sub-leading corrections in the all orders resummation [11] 
(for a pedagogical review, see for instance [12]). The LO BFKL approach can describe HERA 
structure function in a limited kinematical range, i.e. at not so large and small-,x. The 
main assumption at LO is that the processes described through the so-called hard Pomeron 
are driven by gluon ladder diagrams with infinite rungs of s-channel gluons at asymptotic en- 
ergies. A consistent treatment considering higher order resummations is currently in progress 
and should be available soon. An open issue is whether the available energies in the nowadays 
coUiders actually are asymptotic allowing the use of the complete BFKL series. 

In this work we perform a global data fitting to the HERA structure function F2 using as a 
model for the hard Pomeron the finite sum of gluon ladders [13] . The model is based on the trun- 
cation of the BFKL scries keeping only the first few orders in the strong coupling perturbative 
expansion, where subleading contributions can be absorbed in the adjustable parameters. From 
the phenomenology on hadronic collisions [14], just three orders, ~ [a^ ln(l/a;)]^, are enough 
to describe current accelerator data. The hard Pomeron model should be supplemented by a 
soft piece accounting for the non-perturbative contributions to the process. The description 
therefore turns out similar to the two Pomeron model [5], with the advantage of a complete 
knowledge of the behavior on x and Q^. Concerning the energy dependence, the hard piece 
has the logarithmic growth in contrast with the effective powerlike behavior on the model [5]. 
The original model contains a reduced number of adjustable parameters: the normalization Afp 
and the non-perturbative scale /i^ from the proton impact factor and the parameter xo scaling 
the logarithms on energy. More details on those parameter are addressed in the next section. 
In a previous study at Ref. [15], two distinct choices for the soft Pomeron were analised. The 
resulting two Pomeron model was successful in describing data on structure function F2 and 
its derivatives (slopes on and x) for x < 2.5 • 10~^ and 0.045 < < 1500 GeV^. Here, we 
extend the kinematical range of fitting and the models considered for the soft Pomeron. 

This paper is organized as follows. In the next section, we shortly review the main expres- 
sions for the hard piece given by the summation up to the two-rung ladder contribution. In Sec. 
(3), an overall fit to the recent deep inelastic data is performed based on the hard contribution 
referred above supplemented by the remaining soft Pomeron and non-singlet contributions. The 
role played by the hard and soft contributions are investigated in details. There, we also draw 
up our conclusions. 
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2 The hard contribution: summing gluon ladders 



In this section, wc review the elements needed to compute the structure function using the 
finite sum of gluon ladders in the ep collision, with the photon-proton center of mass energy 
labeled by W. Before that, let us shortly motivate the picture of the finite sum of gluon 
ladders. At finite energies, the LLA and NLLA summation implies that the amplitude is 
represented by a finite sum on terms, where the number of terms increases like In s, rather 
than by the solution of the BFKL integral equation. The interest in to take the firts terms in 
the completete series in the truncation is related to the fact that the energies reached by the 
present accelators are not high enough to accommodate a big number of gluons in the ladder 
rungs that eventually hadronize. In the energy range of HERA there can be only a few real 
gluons produced in any scattering event. On the other hand, in the LO BFKL resummation real 
gluons can be emitted without any cost in energy, while in reality the production of a real gluon 
requires an energy equivalent to one or one and half units of rapidity. This violation of energy 
conservation is probably cured in the NLO BFKL resummation or by relying on consistency 
constraint implementing energy conservation in numerical simulations of LO BFKL evolution. 
The truncation of the whole series could be similar to this consistency condition. Corrobating 
the truncation hypothesis, for example the coefficient weighting the term ~ In^ s turn out to 
be compatible with zero considering even the Tevatron data [14] , in contrast with the expected 
from the complete resummation. 

The proton inclusive structure function, written in terms of the cross sections for the scat- 
tering of transverse or longitudinal polarized photons, reads as [16] 

F2(x,Q2) ^ ^^]^aT{x,Q'') + aL{x,Q^)] , (1) 

47r Clem 

^t,l{^, Q') = (Ip / I? ^ ^idk) JC{x, k, k') %{k') , (2) 

where Q is the color factor for the color singlet exchange and k and k' are the transverse 
momenta of the exchanged reggeized gluons in the f-channel. The virtual photon 

impact factor and is the proton impact factor. The first one is well known in perturbation 

theory at leading order, while the latter is modeled since in the proton vertex there is no hard 
scale to allow pQCD calculations. The kernel IC{x,k, k') contains the dynamics of the process, 
for instance, the BFKL kernel. 

The amplitudes can be calculated order by order: for instance the Born contribution coming 
from the two gluon exchange and the one-rung ladder contribution read as, 

T,L{k) %{k), 
cPk cPk' 




k^ k'^ 



^lAk)}C(k,k')%(k'), 



where Os is considered fixed since we are in the framework of the LO BFKL approach. The 
perturbative kernel IC{k, k') can be calculated order by order in the perturbative expansion [16]. 
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The Pomeron is attached to the off-shell incoming photon through the quark loop diagrams, 
where the Reggeized gluons are attached to the same and to different quarks in the loop. The 
virtual photon impact factor averaged over the transverse polarizations reads as [17, 18], 

^'^^ ' 2J0 271 Jo 27T k^pp' + Q^prr' ' 

where p, r are the Sudakov variables associated with the momenta in the photon vertex and 
t' = (1 -r) and p' = (1 -p). 

Gauge invariance requires the proton impact factor vanishing at k' going to zero and is 
modeled in a simple way, 

where A/J, is the unknown normalization of the proton impact factor and fi^ is a scale which is 
typical of the non-perturbative dynamics. These scales will be considered adjustable param- 
eters in the analysis. Considering the electroproduction process, summing the first orders in 
perturbation theory we can write the expression for the inclusive structure function. 



1 /3 a 

+ 



2 V TT X J 



2 

(5) 



In^f 



where the functions I^"'\Q'^, p?) correspond to the n-rung gluon ladder contribution. The quan- 
tity xq gives the scale normalizing the logarithms on energy for the LLA BFKL approach, which 
is arbitrary and enters as an additional parameter. The contributions are written explicitly as. 



/<»' ^ Wm^LjJ'±\ + ^^, (6) 






1 / ^5 \ 

7(2) = 



/131 A , Q'\ 1396 TT* 14,, , 

+ (^+4C(3))ln(^^j+ — -- + -C(3), (8) 

where C(3) ~ 1.202. It is worth mentioning the clear behavior on x and virtuality. The main 
result in Ref. [15] is in a good agreement, in terms of a /doi test, for the inclusive structure 
function in the range 0.045 < < 1500 GeV^, once considering the restricted kinematical 
constraint x < 0.025. The non-perturbative contribution (from the soft dynamics), initially 
considered as a background, was found to be not negligible. In the next section, we perform a 
global analysis, extending the range on x fitted by adding the non-singlet contribution modeled 
through the usual Regge parameterizations and analyzing also the multipole models for the 
soft Pomeron. 
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3 Fitting results and conclusions 



In order to perform the fitting procedure, for the hard piece one uses Eq. (5) and for the soft 
piece we have selected some typical models, as addressed below. First, one takes into account 
the latest version [19] of the CKMT model [20], having a economical number of parameters. 



where A((5^) is the Pomeron intercept. The non-singlet term takes the following form. 



The CKMT model was originally constructed to interpolate between the soft hadronic 
Pomeron phenomenology, where ctpom — 1-08 and the growth for the low proton structure 
function in deep inelastic scattering, where ctpom — 1-2. This is obtained through absorb- 
tive corrections to the bare soft Pomeron, leading to a Q^-dependent intercept. Therefore, in 
phenomenological application considering a two-Pomeron approach, we should be careful in 
verifying the soft Pomeron intercept coming out. 

Another possibility is to select a soft piece which corresponds to the Pomeron with intercept 
equal to one and has the form of non-perturbative truncated In series (soft multipole 
Pomeron) [3, 4], given by 



for the dipole Pomeron and 



A- 



(13) 



Ff\x,Q^') = Q' 



for the tripole Pomeron. 



A ■ 



+ B • In^ 



(14) 



The multipole means the soft Pomeron is considered as multiple poles instead of just a 
single pole. This approach has been successful in hadron-hadron and lepton-hadron (mostly 
vector meson photoproduction) scattering, where the increasing of the cross section on energy 
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Table 1: Parameters of the fit. Procedure (I): hard plus soft terms; Procedure (II): similar 
to analysis (I), but restricting soft Pomeron intercept; Procedure (III): only hard term and 
non-singlet contribution. Procedure (IV): hard -|- soft dipole pomeron terms; Procedure (V): 
hard -|- soft tripole pomeron terms 

can be produced with a unit Pomeron intercept. The multipole model has shown to be quite 
stable in fitting simultaneously hadronic and lepton-hadron data. 

Concerning the hard piece, Eq. (5), we selected the overall normalization factor as a free 
parameter, defined as A/" = |^ Z^e^A/^, considering four active flavors. That contribution 
was also supplemented by a factor (1 — xY" accounting for the large x effects. The remaining 
parameters arc the scales ji^ and Xq as well as the coupling constant (considered at a fixed 
scale in the BFKL formalism). As we will sec in the following, its value presents small variation 
in the fitting, suggesting that it can be considered fixed as ag = 0.2, which is typical in the 
HERA kinematical region. 

For the fitting procedure we consider the data set containing all available HERA data for 
the proton structure function F2 [21]-[30], updated with new HERA experiment [31] and adding 
only the latest (E665 and NMC) data set of fixed target experiments [32, 33]. This choice for 
the selection of data sets follows the procedures used in Ref . [3] . The whole data set contains 
a total number of 1059 experimental points, covering the complete available x and virtualities 
0.045 < < 30000 GeV^. It should be stressed that the inclusion of the older fixed target 
experimental data (SLAC [34] and BCDMS [35]) requires a more sophisticated study for the 
background, for example, following the ALLM parameterization or as in Ref. [3] . However, for 
that purposes the number of free parameters would be significantly larger. 

Here, we have considered the following distinct fitting procedures: 
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Table 2: Parameters of the fit for restricted domain x < 0.07. Procedure (A): hard plus soft 
terms; Procedure (B): hard + soft dipole pomeron terms; Procedure (C): hard + soft tripole 
pomeron terms. 

• (I) Overall fit considering the hard piece and the CKMT soft Pomeron plus reggeon piece, 
Eqs. (9) and (11); 

• (II) Overall fit restricting the soft Pomeron intercept, modifying the expression for A^Q"^), 
see Eq. (10), in the CKMT soft Pomeron. It reads now as, 

• (III) The fit is performed considering only the hard piece plus the non-singlet contribution; 

• (IV and V) Overall fit using the hard Pomeron piece and the soft multipoles Pomeron, 
Eqs. (13) and (14). The main feature is the soft pomeron having an intercept equal one. 

The best fit parameters for these procedures are presented in Table (1). The kinematical 
range covered is < x < 1 and 0.045 <Q^< 30000. In the following we discuss each procedure, 
pointing out the relative contribution from the hard and soft pieces and its quality. 

Lets start by the procedure (I). The parameters for the hard piece remains conssistent with 
the previous analysis in [15], i.e. ctg ~ 0.2, xq c=l 10"^ and ji"^ of order 1 GeV^. Concerning the 

soft Pomeron, the bare intercept Aq coming out quite small in the whole range of Q^. This fact 
seems corroborate a soft Pomeron having an intercept close to one in this two-Pomeron analysis. 
The same happens for the procedure (II), with a modified form for the CKMT intercept. For 
low the intercept comes out small, but reaches at higher virtualities to the upper limit 
A = 0.08. 
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The procedures (IV and V), considering the soft muhipole Pomeron, provide the best quahty 
of fit x^/dof ~ 1. The tripole Pomeron is preferred as a soft background, giving a slightly better 
result than the dipole pomeron. Concerning the hard piece, the value for the fixed coupling 
constant is quite stable, ag 0.2, consistent with the overall HERA value. 

The situation differs only in procedure (III), where a small value for is found, suggesting 
in this case that the best choice would be — 0.119. which coincides with as{Mz). The 
high x^/dof for this case is not quite sizeable, in view of the smaller number of parameters 
considered (8 against 13 from the remaining analysis). 

Let start by the procedure (I). The parameters for the hard piece remains consistent with the 
previous analysis in [15], i.e. ag ~ 0.2, xq ~ 10~^ and /j,^ of order 1 GeV^. Concerning the soft 
Pomeron, the bare intercept Aq comes out quite small in the whole range of Q^. This fact seems 
to corroborate a soft Pomeron having an intercept close to one in this two-Pomeron analysis. 
The same happens for the procedure (II), with a modified form for the CKMT intercept. For 
low the intercept comes out small, but reaches at higher virtualities to the upper limit 
A = 0.08. 

The procedures (IV and V), considering the soft multipole Pomeron, provide the best quahty 
of fit x^/dof ~ 1. The tripole Pomeron is preferred as a soft background, giving a slightly better 
result than the dipole pomeron. Concerning the hard piece, the value for the fixed coupling 
constant is quite stable, ag — 0.2, consistent with the overall HERA value. 

The situation differs only in procedure (III), where a small value for is found, suggesting 
in this case that the best choice would be Og = 0.119, which coincides with as{Mz). The 
high x^/dof for this case is not quite sizeable, in view of the smaller number of parameters 
considered (8 against 13 from the remaining analysis). 

In Figs. 1 and 2, we present the fit result for the inclusive structure function for small 
and large virtuahties. The plots for the different procedures lie on top of each other. In Fig. 
3, we present the relative role of the distinct contributions for the fit. The hard and soft 
pieces are presented separately. We present them explicitly for the virtualities 15, 150 and 1500 
GeV^, where the contributions from fitting (I), (II) and (V) are shown. It is verified that the 
region where the hard Pomeron starts to dominate depends on the vituality. For instance, at 
Q'^ — 15 GeV^ it stays on x ~ 10~^ for the different procedures, whereas is shifted to x ~ 10"^ 
at = 1500 GeV^. A two-Pomeron picture is supported, comparable with those using the 
two-Pomeron analysis in Refs. [3, 5]. 

Additionally, we performed the fit restricting the domain of the Bjorken variable tox < 0.07. 
This procedure is similar to that considered in analysis of Ref . [3] . For this purpose we restricted 
the number of free parameters (equal to 10) for Soft+Hard Pomeron model, as in the Donnachie- 
Landshoff model reanalysis considered in [3]. To this aim, we fixed the constant ag = 0.2 as 
well as the powers Ug = 7 and n„s = 3. The following procedures were considered: (A): hard 
plus soft terms; (B): hard -|- soft dipole pomeron terms; (C): hard -|- soft tripole pomeron terms. 
Results of the fit are presented in Table 2. The quahty of our approach is similar those ones 
found in Ref. [3], at least in this specific domain, x < 0.07. The parameters for the multipole 
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Pomeron seem stable in this kinematical region. There is a change in the case (A), where 
the bare Pomeron intercept is higher than in the procedure (I) producing an usual soft result 
A ~ 0.08. The parameters for the hard piece remain stable in all fitting procedures. 

In conclusion, wc verify that the fitting procedure is equivalent to the model using the 
two- Pomeron approach [3, 5], with the advantage of clear understanding of the behaviors on x 
and of the corresponding hard content. In particular, the behavior on of the hard piece 
sheds light about on further implementations of the hard Pomeron residue in two-Pomeron 
fitting. The fitting procedure shows that the soft Pomeron has an intercept small than the 
usual ajp = 1.08, suggesting that a suitable choice is the multipolc pomeron having intercept 
equal to one. In particular, the tripole pomeron presented the best fit result in all kinematical 
ranges considered here, followed by the dipole one. The fitting using only the hard piece and 
the non-singlet contributions is not completely ruled out, producing a not so high x^/dof result. 
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Figure 1: The inclusive structure function at small virtualities. The procedures (I) - (V) 
produce nearly equivalent curves. The SLAC [34] and BCDMS [35] data are not included into 
the fit. The virtualities are in units of GeV^. 
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Figure 2: The results for the inclusive structure function at medium and high Q'^ virtualities 
(in units of GeV^) 
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Figure 3: The results for the inchisive structure function at = 15, 150 and 1500 GeV^ 
virtuahties with the contribution of soft and hard pieces for the different procedures. 
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